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Abstract 

We discuss the prospects of studying lepton number violating processes in order to identify 
Majorana neutrinos from low scale seesaw mechanisms at lepton-proton colliders. In particular, 
we consider the scenarios of colliding electrons with LHC energy protons and, motivated by the 
efforts towards the construction of a muon collider, the prospects of muon-proton collisions. We 
find that present constraints on the mixing of the Majorana neutrinos still allow for a detectable 
signal at these kind of facilities given the smallness of the Standard Model background. We discuss 
possible cuts in order to further increase the signal over background ratio and the prospects of 
reconstructing the neutrino mass from the kinematics of the final state particles. 
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I. INTRODUCTION 



The first light of the Large Hadron Collider (LHC) has marked the beginning of a new 
era of high-energy collider physics. This actualizes the concept of not only proton-proton 
collisions, but also what kind of physics could be probed in other high-energy collisions 
exploiting proton energies similar to the LHC ones. At the same time, it is the hope that 
new colliders will not only probe the mechanism behind electroweak symmetry breaking 
and the stabilization of the electroweak scale, but also trace and identify the source of other 
known shortcomings of the Standard Model (SM), such as the existence of neutrino masses 
and dark matter. In the case of neutrino masses, the most popular scenario is generally 
considered to be the type-I seesaw models [Mi^ and it has aheady been shown that low-scale 
realizations js ill of such models can give observable signals in collider experiments 12-24|. 

One of the most striking predictions of the type-I seesaw is that neutrinos are Majorana 
particles and thus, that lepton number is broken through their Majorana mass terms. One of 
the most promising probes of the Majorana nature of neutrinos is the search for neutrinoless 
double beta decay (Oi^/3/3), which would produce a distinct peak at the end of the double 
beta decay spectrum (for a recent review see Ref. 25|). Any search for processes violating 
lepton number will typically have a signal suppressed by the tiny neutrino masses 
While in low scale seesaw mechanisms this smallness is due to a cancellation of different 
large contributions, this very same cancellation will suppress any lepton number violating 
signal unless the energy probed in the process is such that some particular contributions are 
enhanced or suppressed due to kinematics. In Oz//3/3 searches, contributions from neutrinos 
with masses above the nuclear scale of ~ 100 MeV are suppressed with respect to those below 
for a discussion) and the cancellation can be avoided. Here we will discuss 



(see e.g. Ref. 
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the prospects of observing the manifestly lepton number violating processes pi^ — )■ i'^jjj- 
Similar to Oi^/3/3, the contributions of mass eigenstates above the energy probed by the 
colliders will be suppressed, evading the cancellation behind the small neutrino masses. 
Moreover, collider experiments, and proton-lepton colliders in particular can, unlike Oz//3/3 
searches, produce these neutrinos on-shell via s-channel contributions with the consequent 
enhancement of the signal. 

The rest of this paper is organized as follows: in Sec. [Tll we review the low-scale seesaw 
model and summarize the current bounds. We then go on to discuss in Sec. IIIIl within the 
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low-scale seesaw model, the lepton number violating processes that constitute our signal. 
Section |IV] is dedicated to providing the details of our numerical computations, while in 
Sec. El we present the results. Finally, in Sec. |Vll we summarize our results and give our 
conclusions. 



II. THE LOW-SCALE SEESAW MODEL AND PRESENT CONSTRAINTS 



Let us consider the Lagrangian of the standard type-I seesaw model, which consists of 
the Standard Model (SM) Lagrangian plus an extra piece containing the allowed couplings 
between the SM fields and additional gauge singlet fermions {i.e., right-handed neutrinos) 



SM 



(1) 



Here, ^ denotes the SM Higgs field, which breaks the electroweak (EW) symmetry after 
acquiring its vacuum expectation value, i'ew- In this work, we will use the basis in which M 
is diagonal with real positive entries. After electroweak symmetry breaking, this produces 
a 6 X 6 mass matrix 

>.^f"'"^V (2) 

where raD = vew'^n, in the basis [i^g, {NkY]'^- In terms of the light mass eigenstates i^* and 
the heavy eigenstates A^*, the neutrino flavor eigenstates can be written as 
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I ~ {5^p - -eaje%)Upiv' + e^iN\ (3) 

where 9 = m^^M~^. In the cases where the heavy neutrinos can be integrated out, this 
implies a non-unitary mixing matrix A/" = (1 — ri/2)U, where rj = 66^ , for the light neutri- 



nos 



29N32| . Inserting this relation into the electroweak interaction Lagrangian, we obtain 
the following coupling between the heavy neutrinos and the charged leptons via the W 
bosons 

^int = -j=W/frOa^N' + H.C. (4) 

The unitarity deviations implied by the heavy neutrinos can thus be used to constrain 
the elements of 9 through the relation r] = 99^ which implies 

l^mT < |r/aa| = 0.005±0.005 (5) 
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at the 90 % Confidence Level (CL). This constraint stems from the effects that such terms 
would have in universality tests of the weak interactions as well as the invisble width of 
the Z [3^. Notice that somewhat stronger bounds \r\f,e\ < 4.0 ■ 10"'^, |?7^^| < 1.6 ■ 10~^ and 
|?7t-,-| < 5.3 • 10~^ are obtained when measurements of Gp from muon decay and the CKM 
unitarity are assumed 3J] instead of the invisible width of the Z. The combination of both 
data sets would, however, result in a larger allowed region since their respective best fits do 
not coincide. Indeed, while the invisible width of the Z prefers non-zero values of 77 in order 
to accommodate the present 2a deviation from the SM, such a result is not favored by the 
measuremets of Gp. For definiteness we will fix here the size of the signal by accommodating 
the data on the invisble width of the Z, setting \r]aa\ = 5.0 ■ 10"'^. However, the rescaling 
of the signal under study to smaller mixing angles is trivial through an overall quadratic 
dependence, as we will discuss later. 

If the masses of the heavy neutrinos display a moderate hierarchy, the signal will mainly 
be dominated by the contribution of the mass eigenstate within the reach of the collider 
energy. Thus, in order to simplify the discussion, we will assume there is only one heavy 
neutrino contributing to the signal, while the other neutrinos, necessary for the generation 
of the observed pattern of masses and mixings at low energies, are heavier and do not 
contribute significantly. The mixing of this neutrino with the SM fermions will be set to 
|6'q,iP = |6'q,P = 0.005 as discussed above. 

Note that, if we also impose the condition that neutrino masses are small while main- 
arg e Yukawa couplings, then an additional constraint |?7q,q,| ■ |?7/3^| ^ |'7a/3p ap- 



tainmg 
plies [S, 
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35[. In particular, strong constraints on the product will then stem from pro- 



cesses such as — >■ 67. Therefore, we will consider scenarios with only one 6^ ^ 0, where 
the mixing is to only one flavour. Thus, we will not consider lepton flavour violation in the 
signal but only lepton number violation. 

Finally, it should be noted that Oz//3/3 decay searches in general provide strong bounds 
on 6ei 271]. However, it is possible to satisfy both the constraints from Oz//3/3 and bounds on 
neutrino masses while still maintaining large Yukawa couplings 11|. 
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FIG. 1: The Feynman diagrams contributing to the process — )• i^jjj at tree level in the 
low-scale seesaw. The left diagram includes a right-handed neutrino in an s-channel, while the 
right diagram has the right-handed neutrino in a t-channel. The left diagram is enhanced when 
the right-handed neutrino can be produced on shell. 

III. LEPTON NUMBER VIOLATING PROCESSES 

Here we will consider the process pi" -^■ to investigate the possibility of producing 

the heavy right-handed neutrinos of the type-I Seesaw mechanism and observe the lepton- 
number-violating signals associated to their Majorana nature at colliders. The advantage 
of lepton-proton collisions with respect to purely hadron or lepton colliders is the cleanness 
of the signal. Indeed, since the lepton number in the initial state is one, the observation 
of a charged antilepton in the final state is a clear signal of violation of lepton number 
by two units, as long as there is no missing energy. This last requisite translates in the 
presence of three jets in the final state, two of them reconstructing an invariant mass equal 
to the W mass. Conversely, lepton number violation in hadronic or colliders implies 
tagging two leptons of the same charge in the final state together with either the presence of 
missing energy in the form of neutrinos or a higher number of jets, making the signal more 
challenging to search for. 

In the type-I seesaw, two diagrams contribute to the signal under study (see Fig. [1]). The 
diagram with the Majorana neutrino in s-channel will be enhanced and dominate the signal 
if the collider energy is high enough to produce it on shell. We will therefore explore this 
possibility, which can help overcome the expected suppression of the signal by the smallness 
of neutrino masses. Thus, the process basically corresponds to the on shell production of the 
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heavy neutrino, via the exchange of a TV boson, and its subsequent decay to i^W~ with the 
W decaying hadronically. The branching ratio for the first decay is roughly 25 % {sgI. On 
the other hand, the production process is weighted by the mixing |6'cp of the heavy neutrino 

with the charged lepton i^- Thus, it is trivial to extend the analysis from the assumed 
values of the mixing angles of |6'ap = 5 ■ 10^'^ to smaller mixings by a general rescaling. 

The electron-proton process we discuss here was already studied in p, 11, [sg], [s^], mainly 
focusing on the DESY experiment but also extended to a combination of LEP and LHC. The 
production cross section and the decay rate of the heavy neutrino were computed in order 
to estimate the rate of the signal that could be expected in a electron-proton collider. In 
Ref. ll| a first numerical simulation of the signal was also performed. The analysis presented 
here is motivated by the improved simulation techniques presently available, in particular 
regarding hadronization processes, and better knowledge of the parton distribution functions 
(PDFs) of the proton. We also extend their analysis to the higher collider energies discussed 



in the LHCeC proposal [38|], which results in an enhancement of the signal and sensitivity to 
higher neutrino masses. We have therefore extended the original analysis to a wider range of 
neutrino masses than initially considered and showed explicitly the dependence of the signal 
on the mediating neutrino mass. 

Furthermore, we also explore for the first time more sophisticated collider technologies, 
such as muon-proton colliders, inspired by the recent efforts towards a muon collider (see 



e.g. Ref. 39[ for a recent status review). This setup would allow to explore different matrix 
elements with respect to the electron-proton option, providing a completely independent 
search channel also with respect to Ou/SP searches. Moreover, the energies that could be 
expected from such a facility imply a huge gain in sensitivity allowing to explore much 
smaller mixings or larger sterile neutrino masses than their electron counterparts. In both 
scenarios we will assume a detector setup based on the capabilities of ATLAS and CMS. 

Regarding possible sources of background, the manifest lepton number violating nature 
of the signal in absence of missing energy makes it very clean and difficult to mimic by SM 
processes. Several sources of background were already discussed in Ref. [ll| where it was 
found that, for the lepton number violating signal under study, the dominant background 
stems from W production with its subsequent decay into the required for the signal^. 



^ Other sources of background discussed in Ref. ll| such as boson-gluon fusion and DIS are negligible 
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The original i~ can be missed for example if it becomes a neutrino or is lost in the beamline. 
We have simulated and studied these two sources of background. In the second case the 
contribution to W production is dominated by the exchange of an almost real photon with 
a very coUinear outgoing electron. In order to compute this process, we have simulated a 
proton-photon collision and convoluted it with the photon PDF in a charged lepton. This 
turns out to be the dominant source of background, in particular when the two final state 



jets are produced by QCD processes. In Ref. 



11| , this background was reduced by means of 



a cut in the invariant mass of the two jets in the final state m^j compatible with that of a ly 
boson Miy, complemented with a cut in the minimum transverse momentum of the outgoing 
{pT,i+)- However, as we will show in Sect. |Vl the presence of neutrinos in the final state, 
necessary to produce the final state positron through SM processes, suggests also the use of 
cuts in missing transverse energy {Ex^miss) to suppress the background. We have found that 
cuts in maximum -Er,miss and minimum pT,e+ actually provide a better signal/background 
ratio than those in rrijj. Moreover, in the case of the p£~ — )■ i'^jjji^ii'i background, two of 
the jets typically originate from a W boson decay. Thus, the invariant mass cut around Mw 
will not avail to reduce this source of background. On the other hand, the cuts in maximum 
ET,miss and minimum Pt/+ prove to be useful for this other background too. This extra cut 

in Ex miss 

thus represents a very useful complementary tool to the cut in pT,e+ in order to 
reduce all sources of SM background. 

IV. NUMERICAL ANALYSIS 

We base our numerical study of the signal and background on the use of the Mad- 



Graph/Mad Analysis 40| software tools supplemented by Pythia 4l| and PGS |42|] to pro 



cess the resulting events. The type-I seesaw description for MadGraph/MadAnalysis was 



obtained via the FeynRules 43| software. For the protons we assume an LHC-like beam 
with an energy of 7 TeV, while several different choices are examined for the lepton beam. 
In the case of an electron beam we study both the conservative setup with a beam energy 



of 50 GeV and a more optimistic setup where the beam energy is 150 GeV 



38 



For muons. 



on the other hand, we consider the beam energies that have been discussed for a future 



because of the lepton number violating nature of the signal considered. 
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Variable Jets Leptons Photons 

PT > 20 GeV > 10 GeV > 10 GeV 

f] < 5 < 2.5 < 2.5 

TABLE I: Acceptance cuts used in MadGraph for our simulations. 



muon collider. These vary between 500 GeV for the more conservative case and 2 TeV for 
the more optimistic proposals 39] . In all cases under study, we employ a detector with 
capabilities similar to the ATLAS and CMS detectors. When generating events, we impose 
the acceptance cuts specified in Tab. [H 

For the simulation of the SM background process pi~ — >■ i'^jjji'i^, we use the default SM 
implementation provided by the MadGraph distribution. In the case of the pi~ — )■ 
which is dominated by the exchange of an almost real photon with a very coUinear outgoing 
electron, we simulate instead the process p7 — ^^jjjv convoluted with the PDF of a photon 
inside the charged lepton. 

For each experimental setup, we simulate the signal for the case of a heavy neutrino 
mass of Mtv = 250, 500, and 750 GeV while setting the mixing to the value required to 
accommodate the invisible width of the Z: |6'q,P = 0.005. This provides a quite optimistic 
signal, but the scaling of the results is trivial with \9a\^, as previously discussed. For all 
cases, we generate a total of 8 ■ 10^ events satisfying the acceptance cuts. At a later stage, 
after presenting the bare results of the simulation, we will implement additional cuts to 
reduce the background. The implementation of these cuts will be made at the MadAnalysis 
level. 



V. RESULTS 



As we are discussing new physics searches, it is fundamental to distinguish between the 
signal and the SM background. Thus, in order to increase the sensitivity of these searches, 
cuts which can suppress the background while not affecting the signal too adversely should 
be investigated. Naturally, an expected background lower than one event would not be very 
worrisome and the constraining factor would then be one of the signal cross section. However, 
as the achievable luminosities for these speculative facilities are uncertain, we will present 
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cuts designed to enhance the signal-to-background ratio regardless of the smallness of the 
cross section. For definiteness, we will assume a baseline integrated luminosity of ~ 100 fb~^ 
and include this value in the plots for reference. Indeed, this value is in the ballpark of the 
10^^ cm~^s~^ being discussed for the LHCeC proposal 38|] for its lower energy version. For 
a 150 GeV electron beam the luminosity could be around an order of magnitude smaller. As 
for a prospective muon-proton collider the luminosities under discussion for a muon collider 
range between 10^^ cm~^s~^ and 10^^ cm~^s~^ 44 1. 

The major question is thus what separates the signal events from the background ones. A 
cut that has been previously studied in the literature is to reject events where the outgoing 

does not have a minimum transverse momentum. The reason for this cut to reject the 
background more strongly than the signal is that, while the signal is mainly due to on-shell 
heavy neutrinos, the background kinematics are such that the positron is predominantly 
emitted in the beam direction. Since all of the SM backgrounds include at least one neutrino, 
which will constitute missing energy, while the signal only includes particles visible to the 
detector, another viable cut could be to impose an upper bound on the missing transverse 
energy of the accepted events. 

In Fig. [21 we show our numerical results for the signal and background distributions in 
both of these variables at an electron beam energy of 50 GeV in order to illustrate how the 
signal and background differ. The behavior for different beam energies, and for the muon 
beams, is similar. In order to illustrate the impact a given cut would have, we show in Fig. [3] 
how the signal and backgrounds for an electron beam would be affected as a function of 
the value at which the cut is implemented. As can be seen from these figures, both cuts 
are indeed effective in increasing the signal-to-background ratio, although there are some 
qualitative differences between them. In particular, due to the peaked nature of the signal 
in the transverse momentum of the outgoing lepton, the cut on this variable is starting to 
deteriorate also the signal quite early in the cases of low Mjy. On the other hand, the shape 
of the transverse energy distribution is similar for all Mn in the area close to zero missing 
transverse energy, thus leading to a uniform behavior for the signal suppression. It should 
also be noted that even a very mild cut on the missing transverse energy would be enough 
to reduce the background by a factor of a few without noticeably affecting the signal and 
thus, at least a modest cut should be implemented in any analysis. For illustration, we will 
implement cuts of Pt/+ > 70 GeV and -ET.miss < 10 GeV. This roughly corresponds to the 
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FIG. 2: Results for the differential cross sections for the 50 GeV electron beam facility. Here, Et is 
the missing transverse energy and px e+ is the transverse momentum of the positron. Blue, red, and 
green lines lines correspond to the signal simulated for heavy neutrino masses of M^r = 250, 500 and 
750 GeV, respectively. The shaded area shows the main backgrounds, the smaller pe~ — t- e^jjjVf^Ve 
in light gray and the larger — )■ e^jjjUf. in green. 



regions where the signals start to be noticeably reduced. 

The impact of the different cuts can be seen from the remaining total cross sections, which 
are presented in Tab. [Tll In the table we also show the cut in the jets invariant mass 



around the mass of the W performed in Ref. 



'33 



ll| for comparison. It can be seen that this cut 



strongly reduces the signal while similar backgrounds suppressions can be achieved through 
the cuts in Pt/+ and -Et.hiIss instead. As can be seen from the table the two cuts are fairly 
independent and complementary in nature, since their combination provides a significant 
background reduction with respect to the implementation of only one of them. Moreover, 
depending on the facility, luminosity and part of the parameter space that is being explored, 
some combination of the two would provide an optimal cut. 

Another interesting kinematic variable is the reconstructed invariant mass of two jets and 
the outgoing lepton Mjj^. In the case of the signal, this should provide a measurement of 
the mass of the on-shell intermediate heavy neutrino, and thus, the majority of all signal 
events cluster around this value. A bump of events at a certain value would not only be 
in favor of signal over background, but also indicate the actual mass of the intermediate 
particle. Even if a search of this sort is performed, the cuts previously discussed still provide 
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FIG. 3: Total cross section for the 50 GeV (upper) and 150 GeV (lower) electron beam facilities as 
a function of a cut on the maximum missing transverse energy (left) and the minimum transverse 
momentum of the positron (right). The blue, red, and green lines lines correspond to simulated 
heavy neutrino masses of M^r = 250, 500 and 750 GeV, respectively. The shaded area shows the 
main backgrounds, the smaller pe~ — ?• e'^jjji^e'^e in light gray and the larger p'j — ?• e'^ jjjve in 
green. The horizontal line represents a cross section of 10~^ fb, which is the required cross section 
to have one expected event at an integrated luminosity of 100 fb~^. 

a complementary suppression of the background. This is illustrated in Fig. HI where we show 
the differential cross section with respect to Mjji at the 50 GeV electron-proton collider as 
an example. An additional cut (not discussed in detail here) which could in principle be 
exploited to enhance the signal over possible backgrounds is the isolation of the final state 
antilepton. 

Finally, in Fig. [S] we show an estimate of the achievable 90 % CL sensitivity to the mixing 
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Cuts M„ [GeV] 


ep [50 GeV] 


ep [150 GeV] 


HP [500 GeV] 


HP [2 TeV] 


No cuts 250 
500 
750 
BG2z^ 
BGlu 


23.53 ± 0.04 
3.351 ± 0.006 
0.1924 ±0.0010 
0.1884 ±0.0004 
5.576 ±0.012 


57.53 ± 0.09 
20.81 ± 0.04 
6.979 ±0.020 
0.9478 ± 0.0020 
27.13 ±0.08 


112.3 ±0.2 
56.32 ± 0.11 
30.35 ±0.07 
3.554 ±0.014 

100.1 ±0.3 


131.9 ±0.2 
113.1 ± 0.2 
76.49 ±0.19 
11.08 ±0.032 
284.8 ±2.2 


PT,e+ > 70 GeV 250 
500 
750 
BG2z^ 
BGlu 


16.65 ±0.05 
2.994 ± 0.007 
0.1721 ±0.0009 
0.0142 ± 0.0002 
1.227 ±0.012 


38.20 ±0.12 
18.6 ± 0.04 
6.272 ±0.020 
0.1381 ±0.0013 
6.65 ±0.04 


80.1 ±0.5 
51.41 ± 0.11 
28.33 ±0.07 
0.8003 ± 0.0061 
29.07 ±0.16 


73.20 ±0.26 
99.49 ± 0.19 
70.16 ±0.20 
3.366 ±0.021 
95.1 ±3.5 


ET,m\Es < 10 GeV 250 
500 
750 
BG2z^ 
BGlz^ 


16.71 ±0.05 
2.182 ± 0.008 
0.1096 ± 0.0007 
0.0045 ± 0.0001 
0.180 ±0.005 


37.48 ±0.11 
11.07 ±0.04 

3.177 ±0.016 
0.0140 ± 0.0004 

0.981 ±0.019 


77.1 ±0.5 
32.39 ± 0.11 
15.09 ±0.06 
0.0373 ± 0.0013 
3.33 ±0.06 


88.46 ± 0.26 
72.69 ± 0.21 
43.90 ±0.20 
0.0691 ± 0.0031 
11.1 ± 1.4 


PT,e+ > 70 GeV 250 
ET,miss < 10 GeV 500 
750 
BG2iy 
BGlz^ 


12.29 ±0.04 
2.018 ± 0.008 
0.1019 ±0.0006 
0.0002 ± 0.0001 
0.1006 ±0.0037 


25.62 ±0.10 
10.24 ±0.04 
2.97 ±0.016 
0.0008 ± 0.0001 
0.546 ±0.014 


57.0 ±0.6 
29.92 ± 0.11 
14.36 ± 0.06 
0.0038 ± 0.0004 
2.05 ±0.05 


54.80 ± 0.24 
66.98 ± 0.21 
41.40 ±0.20 
0.0090 ±0.0011 
8.2 ± 1.3 


75 GeV < rrijj 250 
rrijj < 85 GeV 500 
750 
BG2z^ 
BGlz^ 


1.734 ±0.021 
0.995 ± 0.008 
0.0457 ± 0.0004 
0.0128 ± 0.0002 
0.341 ± 0.007 


3.07 ±0.05 
3.838 ± 0.031 
1.377 ±0.011 
0.0475 ± 0.0008 
1.545 ± 0.022 


5.670 ± 0.2460 

r\ /^o 1 n no 

9.68 ± 0.08 
5.07 ±0.05 
0.1354 ±0.0025 
5.01 ±0.08 


4.960 ± 0.089 
20.80 ± 0.15 
12.67 ±0.13 

0.385 ± 0.008 
13.7 ± 1.6 



TABLE II: The total cross sections in units of fb after applied cuts for the different facihties 
discussed in the text. 

between the heavy neutrino and the charged lepton as a function of the luminosity and for 
different values of the neutrino mass. The sensitivity has been defined as the minimum value 
of the mixing angle that can be excluded with at least 50 % probability at the 90 % CL in 
absence of signal. Given the low statistics expected for small mixings, a Poisson distribution 
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FIG. 4: The differential cross section with respect to the reconstructed mass Mjjg at the 50 GeV 
electron beam facility. The results are shown both without (left) and with (right) cuts on PTe+ 
and i?T,miss- Blue, red and green lines correspond to the signals for heavy neutrino masses of 
Mjsf = 250, 500 and 750 GeV, respectively. The shaded area shows the main backgrounds, the 
smaller pe~ — )• e^jjji^e^e in light gray and the larger p"f — )• e'^jjji'e in green. 

was used to compute the sensitivity. The discrete nature of the Poisson distribution is the 
reason of the abrupt jumps in sensitivity depicted. The sensitivity estimate is conservative 
and based purely on a counting experiment without exploiting kinematic distributions be- 
yond the Et and px cuts such as the reconstructed mass Mjj£ of Fig. HI As can be seen from 
the figure, just based on a counting analysis sensitivities more than two orders of magnitude 
better than the current best fit from the invisible width of the Z would be allowed in these 
facilities for luminosities of 100 fb~^. 

VI. SUMMARY AND CONCLUSIONS 

We have discussed the prospects of testing the existence of heavy Majorana neutrinos 
of a low scale type-I seesaw mechanism in lepton-proton collisions. The question of the 
Majorana nature of neutrinos together with the origin of neutrino masses is one of the 
most fundamental issues still unsolved in particle physics. The quest for answers to these 
unknowns is particularly challenging since all probes testing the Majorana character of 
neutrinos typically face a strong suppression of the signal due to the smallness of the active 
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FIG. 5: The 90 % CL sensitivity to the mixing between the heavy neutrino and electrons (upper 
panels) as a function of the luminosity at the 50 GeV (left) and 150 GeV (right) electron beam 
facilities. The lower panels are the corresponding plots for muons at the 500 GeV (left) and 2 TeV 
(right) muon beam facilities. Red, blue and black lines correspond to heavy neutrino masses of 
Mat = 250, 500 and 750 GeV, respectively. The solid/dotted lines are the sensitivities with/without 
cuts. 



neutrino masses. In low scale seesaw mechanisms, however, this smallness originates from a 
cancellation between different large contributions that could, individually, lead to observable 
signals at searches such as Oi^/3/3 decay or collider experiments. Thus, if the mass spectrum 
of the heavy Majorana neutrinos is such that some contributions are enhanced or suppressed 
by the kinematics of the process, the cancellation responsible for the smallness of neutrino 
masses will not take place. This is particularly true in collider searches, as the one discussed 
here, given the fact that the contribution of a given neutrino can be enhanced via its s- 
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channel production if the colhder has enough energy to produce it on-shell. 

We have simulated the lepton number violating signal that the contribution of a heavy 
Majorana neutrino within the reach of the lepton-proton collider would provide. The sim- 
ulations were performed for two types of prospective facilities. We have extended previous 
analyses of electron-proton colliders and updated them to the recent LHCeC proposals in 
which 7 TeV protons collide against electrons with energies between 50 and 150 GeV. As 
a more ambitious facility we also considered a muon-proton collider inspired by the efforts 
towards a muon collider. We simulated muons with energies between 0.5 and 2 TeV colliding 
against 7 TeV protons. Such a facility would not only allow to probe smaller mixings and 
higher neutrino masses through higher cross sections of the signal at higher energies, but 
also provide a complementary search channel, since it would probe areas of the parameter 
space completely independent to the electron-proton collisions or Oz//3/3 searches. 

In order to simulate the signal, we have set the mixing between the new heavy neutrinos 
and the charged leptons to the value required to accommodate the 2cr deviation of the 
invisible width of the Z with respect to the Standard Model prediction. Even this optimistic 
assumption leads to a small signal, as expected from the challenging task of testing Lepton 
number violation. We have shown that, despite the smallness of the expected cross section, 
the facilities under study would allow for an observable signal in areas of the parameter space 
still allowed by present data. Moreover, the signal from these facilities is particularly clean, 
with a very small Standard Model background to obscure it. The dominant background 
is originated from W production, with its subsequent decay into an and the original 
t~ missed by the detector. Nevertheless, we have studied how the signal over background 
ratio could be improved further through different kinematic cuts so as to allow the search 
of even smaller values of the mixing angle, provided sufficiently high collider luminosities. 
A conservative estimate of the sensitivity to the mixing between the heavy neutrinos and 
charged leptons yields an improvement of more than two orders of magnitude with respect 
to the present constraints for luminosities of 100 fb~^. 

We conclude that lepton-proton colliders provide a particularly clean probe of the elusive 
Majorana nature of neutrinos and the type-I seesaw mechanism and would complement 
other ongoing searches by exploring different parts of the allowed parameter space. 
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